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Abstract. We give sufficient conditions for the existence of positive travelling wave 
solutions for multi-dimensional autonomous reaction-diffusion systems with distributed 
delay. To prove the existence of travelling waves, we give an abstract formulation of 
rS^ • the equation for the wave profiles in some suitable Banach spaces, and apply known 

C^ , results about the index of some associated Fredholm operators. After a Liapunov- 

Schmidt reduction, these waves are obtained via the Banach contraction principle, 
as perturbations of a positive heteroclinic solution for the associated system without 
diffusion, whose existence is proven under some requirements. By a careful analysis 
^ ' of the exponential decay of the travelling wave profiles at — cx), their positiveness 

>0 , is deduced. The existence of positive travelling waves is important in terms of 

1^ ' applications to biological models. Our method applies to systems of delayed reaction- 

rvq ■ diffusion equations whose nonlinearities are not required to satisfy a quasi-monotonicity 

ly-T I condition. Applications are given, and include the delayed Fisher-KPP equation. 

o 
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1. Introduction 

For the last decades, there has been an increasing number of studies in travelhng wave 
fronts for delayed diffusion equations, and several methods to prove their existence have 
been developed. 

In this paper, we are concerned with the existence and positiveness of travelling 
waves connecting two equilibria, for a class of A^- dimensional systems of reaction- 
diffusion equations with distributed delay in the reaction terms, of the form 

Ou 

— (t, x) = Au{t, x) + /(mj(-, x)), teR, xe W. (l.l) 

Here, f : C := C([— r, 0];M^) — )> M^ is continuous, C is equipped with the norm 
WfWoo = supgg[_^o] 19^(^)1) fo^ some fixed norm | ■ | in R^, and r > 0. As usual, 
Ut{-,x) denotes the restriction of a solution u{t,x) to the time interval [t — r, t], i.e., 
Ut{0,x) = u{t + 6,x) for —T < 6 < 0,x E MP. For simplicity, we consider all the 
diffusion coefficients equal to 1 in (1.1), but all our results apply to the more general 
case of the diffusion term given by DAu(t, x), where D = diag {di, . . . , dj\f) with di > 0. 

We are mostly interested in situations where (1.1) represents a population dynamics 
model, or another biological model. Typically, we want to obtain conditions for 
the existence of a travelling front connecting two steady-states, zero and a positive 
equilibrium K G M^. Due to the biological interpretation of the model, only non- 
negative solutions are meaningful, therefore we look for positive travelling wave solutions, 
connecting to i^ as t goes from — oo to oo. 

With the method presented here, such positive travelling waves are obtained 
for large wave speeds, as perturbations of a positive heteroclinic solution for the 
corresponding functional differential equation (FDE) without diffusion, 

u'it) = f{ut), t G M (1.2) 

(where Ut E C denotes the function Ut{0) = u{t + 6) for 6 G [— r, 0]), whose existence 
we shall prove under some requirements on /. This idea is not original, and has been 
exploited in the literature (see e.g. [9, 10]). When compared with [9], in the present 
paper the major novelty is that we give conditions for the travelling waves to be positive. 
We also note that [9] considers delayed reaction-diffusion equations with a global space 
interaction, a situation not considered here, for the sake of simplicity. Our results can 
however be extended easily, to take into account non-local effects. On the other hand, 
[10] deals with scalar reaction-diffusion equations with one single discrete delay of the 
form 

du u 

where f{u(t, x),u{t — r, x)) = —u(t, x) + g{u{t — r, x)) and g : [0, oo) — )■ [0, oo) is C^- 
smooth, (7(0) = 0,g{K) = K for some K > Q., and 5''(0) > 1. Assuming that the 
two equilibria and K are hyperbolic, under some further assumptions the existence 
of positive and in general non-monotone travelling wavefronts connecting to -ft' was 
established in [10]. 
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Recently, several techniques have been developed to prove the existence of travelling 
wave fronts for delayed diffusion equations. They are often based on the application of a 
fixed point theorem in an adequate Banach space, which requires a quasi-monotonicity 
condition, either for the original equation (1.3) [16, 23] or, more recently, for some 
auxiliary equations [17]. These methods are usually combined with a monotonic iteration 
scheme, associated with the construction of a pair of upper and lower solutions. See 
[16, 17, 23] and references therein. We emphasize that our method applies to systems 
(1.1) with non-monotone nonlinearities, in the sense that we do not impose on / any 
type of quasi-monotonicity condition, as defined in [20, 23]. 

Before introducing our hypotheses, we set some standard notation. For d = 
{di, . . . ,diy) G M^, we say that d > (respectively d > 0) ii di > (respectively 
di > 0) for i = 1, . . . ,N. In C, we consider the partial order (p > ip ii and only if 
(f){9)-ij{9) > for ^ e [-r, 0]; in a similar way, (p > tp ii (p{9)-ij{9) > ior 9 e [-r, 0]. 
As usual, C+ denotes the positive cone C{[—t, 0]; [0, oo)^). 

For /, the following hypotheses will be considered: 

(HI) /(O) = f{K) = 0, where K is some positive vector; 

(H2) (i) / takes bounded sets of C into bounded sets of R'^ and is C^-smooth; 
furthermore, (ii) for all M > there is /3 > such that /i(v?) + /3v^j(0) > 0, 
i = l,...,N,ior a\\(f eC with < v? < M; 

(H3) for Eq. (1.2), the equilibrium u = K is locally asymptotically stable and globally 
attractive in the set of solutions of (1.2) with initial conditions ip G C+, (/^(O) > 0; 

(H4) for Eq. (1-2), its linearized equation about the equilibrium has a real 
characteristic root Aq > 0, which is simple and dominant (i.e., ^z < Xq for all 
other characteristic roots z); moreover, there is a characteristic eigenvector v > 
associated with Aq. 

We summarize the main results in this paper as follows. In Section 2, we assume 
(H1)-(H4) and establish the existence of a positive heteroclinic solution u*{t) to (1.2), 
with u*{—oo) = 0,u*{oo) = K, and asymptotic behaviour 0{e^°'^) at — oo. In Section 
3, for large wave speeds we prove the existence of travelling wave solutions for (1.1), 
connecting to K. The profiles of these waves are obtained as perturbations of u*{t) 
via a contraction principle argument. For this, we generalize the procedure in [9], 
and give an abstract formulation of the wave profiles as solutions of an operational 
equation, acting in suitable Banach spaces, which incorporate a desirable exponential 
decay 0{e^^) at — cxd, < /x < Aq. Some nice results of Hale and Lin [13] on the index 
of some associated Fredholm operators are used, and a Liapunov-Schmidt reduction 
effected, to set up the right framework for the application of a contraction principle. 
As mentioned above, an existence result of travelling waves connecting two hyperbolic 
equilibria was already obtained in [9], for a class of reaction-diffusion equations with 
global response, but for such waves no exponential decay at — oo was derived in [9], nor 
their positiveness. By a careful analysis of the behaviour of the wave profiles at — oo. 
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in Section 4 we prove that there are positive travelhng waves if the wave speed is large 
enough, and exphcitly give their asymptotic decay at — oo. Section 5 is dedicated to 
apphcations, which include the Fisher-KPP equation with delay and a 2-dimensional 
chemostat model. An important theorem on the asymptotic behaviour of solutions of 
perturbed linear autonomous ordinary FDEs is given in the Appendix. This result 
generalizes to the case of FDEs with distributed delay a result by Mallet-Paret [18] , for 
FDEs with discrete time-delays (or time-shifts), and is often used in Sections 2 and 4. 

2. Positive heteroclinic solution for Eq. (1.2) 

In this section, we prove the existence of a positive solution u*{t) of the ordinary FDE 
(1.2) connecting the equilibrium to the positive equilibrium K. We recall that a 
function u{s) defined on a set S and with values in R^ is said to be positive if all its 
components 'Ui(s), . . . ,un{s) are positive functions on S. 

Theorem 2.1 Assume (H1)-(H4). Then: 

i) There exists a heteroclinic solution u*(t),t G M, of Eq. (1.2), with u*{—oo) = 
0,M*(oo) = K; 

(a) u*(t) is positive, t G IR; 

(Hi) u*{t) = ce'^"*v + 0(e*^^'*'"~^^*) at — oo, for some c > and each fixed e > 0. 

Proof, (i) Consider the linearization of (1.2) about 0, 

u'{t) = Lut, where L = Df{0), (2.1) 

and its characteristic equation 

det Ao(A) = 0, where Ao(A) = L(e^-J) - XL (2.2) 

Recall that A is a solution of (2.2) if and only if A G cr(A), where A is the infinitesimal 
generator associated with the semi- flow of (2.1). 

Let Ao > be the leading (simple) eigenvalue of (2.1) given in (H4), and v G 
IR ,v > 0, such that Ao(Ao)v = 0. Choose 7 > with 7 < Ao < 27 and such that the 
strip 7 < i?e A < Ao does not contain any root of (2.2). Define Xo(^) = e^"^v, 9 G [— r, 0], 
and decompose the phase C as C = PQ)Q, where P =< xo > and Q is the complementary 
space given by the formal adjoint theory of Hale [14]. Then there are neighbourhoods 
A^o, A^i of in P, Q, respectively, and a C^ map w : Nq ^ Ni with w{0) = 0, Dw{0) = 
such that the local 7- unstable manifold of for Eq. (1.2) is given by 

W^(0) = {0 + w(0) :0G A^o}. 

Note that ip G 14^(0) if and only if there is a full trajectory Ut = Ut{ip) {t G M) of (1.2) 
with Mo = V^, Mt G Ao + A^i for t < and u{t)e~"'^ — )■ as t — )■ —00. See Krisztin et al. 
[15], Hale and Lunel [14, Sec. 10.1-10.2], and Diekmann et al. [5, Sec. 8.4]. 

We now argue as in [19]. Write w{t) = {wi(t), . . . , WN(t)),\ = (vi, . . . , vat). Since 
Dw{0) = 0, then lim||0|j_>.o,(/,eAfo yVy-^ = and we have limc_s.o 1?" = 0. Thus, there 
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is Co > such that 1^1(0x0)100 < ce~'^'''^Vj/2 for c G (0,Co],i = 1, . . . ,iV, which imphes 
that for < c < Co we have 

mill fce^»'v, + u;,(cxo)(^)) > ce-^«"v,/2 > 0, i = l,...,N, (2.3) 

and therefore cxo + u]{cxo) ^ W{0) fl C_|_ for all c G (0,Co]. Fix e.g. c = co, denote 
= coXo + 'w{coXo) and consider the full trajectory u^ = Ut{(f)), t E M. We have 
u^ G W{0) for t < 0, hence m^ has the form u^ = c(t)xo + 'w{c(t)xo)- Since the map 
t h^ u^ and the canonical projection of C on P are continuous, c(t) is continuous as well, 
with c(t) — 7- as t — 7> —00. This implies that there is T < such that c(t) < Cq for t < T. 
On the other hand, if c(to) = for some to < T, then u^^ =0, which is not possible. 
From (2.3) it follows that u*{t) > for t < T. Now, from (H3) we have u*{t) -)■ K 
as t — > 00. This means that u*{t) is a heteroclinic solution of (1.2) connecting the two 
equilibria 0,i^, with u*{t) positive on some interval (— oo,T). 

(ii) Choose M > such that u*{t) < M, t e R, i = 1, . . . , N. For the sake of 
contradiction, suppose there is t > T and i G {1, . . . , A^} := I with u*{t) < 0. Define 
t* = min{t > T : u*{t) = for some j G /} and take i G / such that u*{t*) = 0. For 
M as above, let (3 be as in (H2), i.e., fj{(p) + l3ipj{0) > 0, for j e I and < ip < M. 
Writing u*{t) in integral form, 

t 



um= e-^^'-^\f,{u:)) + /3u*{s))ds, te 



we obtain 



0= r e-^(**-^)(/,«)) + /3<(s))rf5. 



where /«(«*) + f3u*{s) > for s < t*. Hence /«(«*) + f3u*{s) = for s < t*, and in 
particular u* satisfies the scalar ODE y' = —(5y for s <t* . Thus u*{s) = for s < t*, 
which is not possible. 

(iii) We note that ul belongs to W{Q) for t < 0, thus -u* = 0{e"'^) at —00 and u*{t) 
satisfies u'{t) = Lut + h{t), with h{t) = f{ul) - Lu^ = 0{e^'^^) at -00. From Theorem 
6.2 (see Appendix), for each e > we deduce that u*{t) = z{t) + 0(e*^^'''~^^*) at —00, 
where z{t) = ce^^^\ for some (positive) c G M. Thus, u*{t) = 0(6^*"'*) at —00. ■ 

Remark 2.2 In fact, one could use [10, Lemma 4] and its constructive proof to derive 
that there is a complete solution u*{t) of (1-2), with u*{—oo) = 0,u*{oo) = K, and 
u*{t) > for t < 0. This proves assertion (i) of Theorem 2.1. In order to prove that 
u*(t) = 0(e'^"*) at —00 it is however more convenient to explicitly construct u*{t) as a 
perturbation of the eigenfunction e'^°*v as above. This asymptotic result will be crucial to 
prove the existence of positive travelling waves for (1.1), if the wave speed is sufficiently 
high. On the other hand, if we assume that the interior of the positive cone C+ is 
positively invariant for the flow of (1.2), as an alternative to hypothesis (H2)(ii), then 
the positiveness of u*{t) on M follows immediately from the fact that u*{t) is positive in 
the vicinity of —00. 



Positive travelling fronts for delayed systems 6 

3. Existence of travelling waves and their asymptotic decay at — oo 

Throughout this section, for simplicity we assume that (H1)-(H4) are fulfilled, but in 
fact some of the hypotheses can be weakened (cf. Remark 3.12). We shall prove the 
existence of travelling waves for (1.1) which will be obtained as perturbations of the 
heteroclinic solution u*{t) of (1.2). The asymptotic behaviour at — oo of u*(t) given 
in Theorem 2.1(iii) will be important to study the asymptotic decay of such waves at 
— oo; however its positiveness is irrelevant here, and will be only used for the analysis 
in Section 4. 

For a unit vector w G W\ we look for wave solutions of (1.1) with direction 
w and speed c > 0, connecting the equilibria to K, i.e., solutions of the form 
u(t, x) = (f){ct + w ■ x) with 0(— oo) = 0, <?i>(oo) = K. 

The equation for the travelling wave profile is given by 

0"(t) - c0'(t) + /,(0,) = 0, teR, (3.1) 

where fc{4>) = /(0(c-)), with subject to the conditions 

0(-oo) = 0, 0(oo) = K. 
With e = 1/c, (3.1) is equivalent to 

eYit) - 0'(t) + /(0,) = 0. (3.2) 

We also consider Eq. (3.2) with e = 0, in which case it reduces to Eq. (1.2). 

Let Cb(M, M^) be the space of all continuous and bounded functions from M to M^, 
with the supremum norm |||/||oo = sup^gjj |l/(s)|. As a particular case of the framework 
in [9], we have the following result: 

Theorem 3.1 [9] Let f have the form f{(f)) = F{(f){0),g{N(f))), (/) E C, for some bounded 
linear operator N : C ^ M.^ and g : M.^ -^ M.^ , F : R^^ -^ M^ C^ -smooth functions. 
Suppose also that: 

(i) /(O) = f{K) = for some K e R^ , 

(a) for Eq. (1-2), the equilibrium u = is hyperbolic and unstable, and the 
equilibrium u = K is locally asymptotically stable. 

Then, if there is a heteroclinic solution u*{t) for (1.2) connecting to K , for each 
unit w eRP there are a neighbourhood V of u*{t) in (^^(IR, M^) and a constant c* > 0, 
such that for c> c* the set of travelling waves for (1.1) in V , with direction w and wave 
speed c, constitutes a C^ -manifold of dimension m, where m is the dimension of the 
unstable space for u(t) = Df{0)ut. 

In this section, the idea is to retrace some arguments in [9] for the proof of Theorem 
3.1 adapted to the case of (1.1), but in appropriate Banach spaces, which will allow us 
to deduce not only the existence of travelling wave solutions for (1.1), but also their 
asymptotic behaviour at — oo. This behaviour will be used in Section 4, to prove the 
existence of positive travelling waves. 

In addition to Cb := Cfe(]R, M^), we introduce the following Banach spaces: 
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R^) = {y e Cb'.y' e Cb} with the norm \\y\\i = ||y||oo + \\y'\ 
Co = {y E Cb '■ \im.s_^±ooy{s) = 0} is considered as a subspace of Cb] 



ci-.-- 



ci( 



Cq = {y E Cl : y,y' G Co} is considered as a subspace of C^; 
C^ = {y E Cb '■ sup^<Q 6^^**11/(5) I < 00} (for /i > 0) with the norm 



\y\ 



max{||i/||oo, IblL} where \\y\ 



supe 

s<0 



-/IS 



y{s) 



\y\ 



\y 



IM' 



Cj^ = {y e Cl ■.y,y'e C^}, with the norm \\y\\i,f, 
C^fl = C^ n Co is considered as a subspace of C^. 
By the change of variables 0(t) = w{t) + u*{t), (3.2) becomes 

£V'(t) - w'{t) - w{t) = -w{t) - Df{ul)wt - G(£, t, w), te 

where 

C(5, t, ti;) = /(t^, + ul) - f{ul) - Df{ul)wt + e^u*"{t), 

subject to the conditions w{—oo) = w{oo) = 0. The roots of the characteristic equation 
associated with e^w"(t) — w'(t) — w{t) = are 

. . 1 - VI + 4^2 



(3.3) 
(3.4) 



a[e 



m 



1 + VI + 4^2 



and satisfy Q;(e) — )■ — l^,/3(e) — )■ 00 as £ — ?■ 0+. In the case of different diffusion 
coefficients di > 0,i = 1,...,N, instead of a{e),f3{e) one has to consider ai{e) , f3i{e) , 
the solutions of the characteristic equations die'^z^ — z — 1 = 0,i = 1, . . . ,N, but the 
arguments are similar (cf. [9]). 



A bounded function w : 



-)> 



pN 



is a solution of (3.3) if and only if 



Jw{t) = H{e,w){t), te 



(3.5) 



where 



and 



VI + 4^2 



Jw{t) = w{t) 



Hie,w)it) 
t 



<'-'\w{s) + Df{u:)ws]ds 



=a{£){t-s) 



= -(t-^) 



. IVTT^ 

»(^)(t-')G{6,s,w)ds + 



{w{s) + Df{u*)ws)ds- 



Jie)(t-s) [^(g) ^ Df{ul)ws + C(e, s, w)] ds 



Our purpose is to apply a contraction principle argument in order to obtain a 
solution of Eq. (3.5), for £ > small and w close to 0, in adequate spaces C^. We first 
analyse the linearity J, by introducing some auxiliary equations and operators. 

Define 

{Ty){t) = y\t) - DfiuDy,, y G C,\t G R. 

We easily see that J : Cq — )■ Co, T : C^ — )■ C^ are linear bounded operators and 
w I— 7- H{w, e) maps Co in Co, for e > (cf. [9]). For yU > 0, we also define 



T • = 

■J-fj, ■ 



TU : a. ^ C, 
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Lemma 3.2 Let fi > 0. Then T^ and J\c „ : C^,o ~^ C^^q are bounded operators. 

Proof Since the map t \-^ \\Df{ul)\\ is continuous for t G IR and \\Df{ul)\\ -)■ ||D/(0)|| 
as t -^ -oo, \\Df{u;)\\ -^ \\Df{K)\\ as t ^ oo, then M := sup^gj^ ||D/(m*)|| < oo. It 
follows that \{Ty){t)\ < max(l, M)||?/||i,^ for t > and e-^*|(Ti/)(t)| < max(l, M)|||/||i,^ 
ioTt<0,yeCl,. 

For y e Cf,, we now have \{Jy){t)\ < (2 + M)||t/||^ for t > and e-^*|(Jt/)(t)| < 
[1 + (1 + M){fi + l)-i] \\y\\^, hence J(C^,o) C C^,o. ■ 

Consider the linear variational equation around the heteroclinic solution u*{t), 

y\t) = Df{ul)yt. (3.6) 

Define the operators L(t) := Dfiul) in (3.6), with L(— oo) = Df{0) and L(oo) = 
Df{K). Hence, Eq. (3.6) is asymptotically autonomous, with limiting equations (2.1) 
and y'(t) = Df{K)yt, respectively at — oo and oo. 

Lemma 3.3 Consider fi G (0, Aq) such that there are no characteristic roots A of (2.1) 
with 5R A = /i. For T^ : Cj^^ C^ defined as above, 

Im{T^) = C^, dim Ker (T^) = r^, 

where r^ = #{A G C : det Ao(A) =0,^\> fi}. In particular, r^ = 1 for fi close to Aq. 
Moreover, Ker (T^) C C^,o- 

Proof Clearly, equation y'(t) = Df{K)yt is asymptotically stable, and the autonomous 
equation (2.1) admits a "shifted exponential dichotomy" in M with the splitting made 
at fi and exponents fi — 6, fi + 6, for S > small. See Hale and Lin [13] for definitions, 
and note that C^ = C°(/i, 0) in the notation in [13]. From [13, Lemma 4.3], there is 
T > such that (3.6) has a shifted exponential dichotomy on (— oo, — T] and [T, oo). 
We now apply Lemma 4.6 of [13] to (3.6). It follows that T^ is a Fredholm operator, 
with index Ind{T^) given by 

Ind{T^) = dim I m{p- {-t)) - dim/m(P+(t)), t > T, 

where P^ {—t) , P^ {—t) and P^ (t) , P^ (t) (t > T) are the projections associated with 
the (shifted) exponential dichotomies for y'{t) = Df{fS)yt and y'{t) = Df{K)yt, 
respectively. From [13, Lemma 4.3], we also have that P^{—t) — )• P^,P^{t) — )■ P+ 
as t — 7- oo, where P~ is the canonical projection from C onto the /i-unstable space 
E~ for y'{t) = Df{0)yt, and P^ is the canonical projection from C onto the unstable 
space E^ for for y'{t) = Df{K)yt. We have E^ = {0} and dimii^~ = r^, where 
r^ is the number of characteristic values for (2.1) (counting multiplicities) with real 
parts greater than fi. Hence Ind{T^) = r^. On the other hand, the index of T^ is 
defined by Ind{Ty) = dim Ker(T^) — codim/m(T^). Again by [13, Lemma 4.6] we have 
dim i^er(T^) = dim.E~ = r^, yielding that Im{T^) = C^. 

For y G Ker{Tn), from the definition of shifted exponential dichotomy we have 
hmi^oo y(t) = 0. Thus, Ker (T^) C C^,o- ■ 
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Similarly to what was done for T, we now restrict the domain and range of the 
operator J. With D := d/dt + id, consider the commutative diagram 

It is easy to check that this diagram is well defined, and that D is one-to-one and 
surjective. Since T^ = D o J is surjective, we may conclude that J is also surjective. 
Moreover, 

Lemma 3.4 Consider ^ G (0, Aq) such that there are no characteristic roots A of (2.1) 
with ^\ = fi. Then, for the operator J\cf,,o '■ C^i,o ~^ C^^ we have Ker {J\ci^o) = 
Ker (T^) and Im (J|c^_o) = C'm,o- 

Proof Recall that Ker (T^) C C^,o- Clearly, for w G C^^ we have Jw = 
if and only if w'(t) = Df{u^)wt, and then w' G C^. We therefore deduce that 
{Ker J)nC^ = (Ker J) n Cj,, and Ker (J|c„„) = Ker (T^). 

We now prove that Im{J\c^o) = C^^. Indeed, for y G C^^ we have that 
C, := y — Jy E Cjj^ and D^{t) = y(t) + Df{ul)yt, hence DC^ G C^,o- Equation Jw = y is 
equivalent to J{w — y) = C,, and therefore it possesses a solution x = u! ~ V ^ Cj^- After 
applying D to both sides of the latter equation, we get T^x = DC, G C^,o- Since the 
cj-limit operator T^(oo) is hyperbolic, we may invoke Lemma 3.3 from [9] to conclude 
that x(oo) = 0. Thus w{oo) = 0, and J : C^^ — )■ C^^ is surjective. ■ 

We now focus our attention on the non-linearity H of Eq. (3.5). Proceeding as in 
[9], one sees that H{w, e) G Cq for each e > and w G Cq. We want however to consider 
the maps if (■, e) restricted to some neighbourhood of zero in C^^, for e > and /i > 0. 
We start with an auxiliary lemma: 

Lemma 3.5 Let X,Y be normed spaces and C G O G X . Suppose that the set C is 
compact, O is open and F : O ^ Y is a continuous map. Then for every a > there 
exists 6 > such that 

\F{x + z)- F{x)\ <a, X gC, \z\ < 6. 

Proof. By the continuity of F, for each x G C there is 6{x) > such that if \z\ < 26{x) 
then X + z G O and \F{x + z) — F{x)\ < a jl. Since C G V}x(^KBi(x){x) is compact, 
there is a finite subcover {-B5(xj)(a;j)}^i of C. For each x G Bs(xj){xj) fl C and \z\ < 
6 ■.= mm{6{xj)}, we have \F{x + z) - F{x)\ < \F{x + z) - F{xj)\ + \F{x) - F{xj)\ < a, 
which proves the lemma. ■ 

Lemma 3.6 Assume (H1)-(H4) and consider fi G (0, Aq). For any 6 > 0, there are 
e* > (independent of fi) and a > such that H{w,e) G C^^ for any e > and 
w G C^,on5^(0), and 

\\H{w,e)\\,<5i\\w\\^ + l), 

\\Hiw,s)-Hiv,e)\\^<5\\w-v\\^, for w,v G C^,or] 3^(0), e G (0,6*) (3.7) 

where B^{0) is the a -neighbourhood of in C^. 
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Proof. We write H = Hi + H2 + H^, where 

gO!(e)(t-s) 
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Hiie,w)it) 
H2{e,w){t) 
H3{e,w){t) 



VI + 4e2 



— e 



its) 



{w{s) + Df{ul)ws)ds, 



Vl + 4e2 



—00 
00 



;/3(£)(t-«) [^(s) + Df{ul)ws + ^(e, s, w)] ds, 



and G is given by (3.4). Let M = sup^gjpj^ ||D/(-u^)|| as before. For t E M., e > and 
/i > 0, we have 



^a{s)(t-s) 



<t-s) 



e^'ds 



-Me)(t-s) _ g-(t-s) 



e^^'ds 



Vl + 4^2 
< , ^ /" e"(^)(*-^)(l - Vl + 4£2) + Vl + 4^2 f e' 

" VI + 4£2 y_^ V 

(Vl + 452) - 1) / e"(")(*-')e^^ ds + Vl + 4^2 I (e"(=)(*-^) - e-^*"^)) e^' ds 
J —00 J —00 



Vl + 4^2 

1 

~ Vl + 4^2 

1 

~ Vl + 4^2 
where 



2Vl + 4e2 - 1 ^ Vl + 4£2 
/i + 1 



;U — a{e) 



oM-t 



Vl + 4^2 - 1 (l + a(£))Vl + 4£2 



e^ < Ci(£)e 



^t 



(3.^ 



Ci(^) 



« £ 



1 + a(e) -> as £ -)■ 0+. 



Vl + 4^2 
From (3.8), we obtain 

\\Hi{e,w) - Hi{e,v)\\^ < Ci{e){l + M)\\w-v\\^, w,veC^,e> 0. 

Since Hi{e, 0) = 0, in particular Hi{e, w) G C^ for w E C^ and e > 0. 
For < /i < /3(e) and t G M, we now have 



(3.9) 



,a(£)(t-s)gM«^g 



ofJ-t 



+00 



,/3(e)(t-s)gM«^g 



sM* 



Consider e.g. M^ equipped with the maximum norm. For t G M, e > 0, w,f 

i = 1, . . . , A^, we have 

\G,{e,t,w)\ < e'\u*"{t)\ + \f,{wt + <) - /,«) - Df,{u;)wt\ 

< e'\u*"it)\ + ||D/,« + ^,,twt) - /^/.«)|| Iktiloo 
and 

\Gi{e,t,w) -Gi{e,t,v)\ < \\Dfi{vt + ul + 9i^t{wt - Vt)) - Dfi{ul)\\\\wt - Vt\\oo, 
for some ^j^t, ^j^t G (0,1) for t G M. 



(3.10) 

(3.11) 
(3.12) 
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Note that u*" G C^ for < /i < Aq. In fact, u* E C^ from Theorem 2.1, 
hence Eq. (1.2) and the smoothness of / lead to |'U*'(t)| < Mo||mJ'||oo5 from which 
we derive Hm*']]^ < Mo||m*||^, for some Mq > 0. By differentiating, we obtain 
u*"{t) = Df{ut){u*%, thus \\u*"\\, < M\\u*'\\,. 

In order to simplify the notation, for each fi,a > write C^^ ni?o-(0) to denote the 
(T- neighbourhood of in C^,o- Since u* is uniformly bounded on M and / transforms 
bounded sets of C into bounded set of R^, then u*' is uniformly bounded on M 
and u* uniformly continuous on R. Thus, /C = {u^, t G R} C C is compact. The 
continuity of Dfi : C — ?> C{C,C) and Lemma 3.5 imply that, for each 6 > fixed, 
there is a = a{6,fi) > such that \\Dfi{vt + u^ + 6i^t{wt — Vf)) — Dfi{ul)\\ < 6 for 
w,v e C^,o n 5^(0), t G IR. From (3.11), (3.12), we get 

\G{e,t,w)\ < e^\u*"{t)\ + 6\\wt\\oo, u^ G C^,o n 5,(0), 
\G{e,t,w)-G{e,t,v)\ < 6\\wt - vt\\oo, w,v e C^,o n fi<,(0). 

From these estimates and (3.10), we conclude that H2{e , w) , H^i^e , w) G C^ for all 
w e Cf,n B^{0) and e > 0, /i G (0, Aq), with 

\\H2{e,w)\\,<- -^-^==={e'\\u*"\\, + S\\w\\,), 

(/i — a(e))vl + 4:6^ 

\\H,ie,w)\\,<—- L^^\e^u*"\\, + il + M + 5)\\w\\l (3.13) 

(/3(e) -;u)v1 + 4^^ L ^ 

Furthermore, for \\w — v\\^,w,v G C^,o H -8,(0), e > 0, we get 



Y) 

l + M + 6 
{/3{e 
On the other hand 



H2{€,W)-H2{€,V)\\^< —-j====\\w-v\\^, 

[fi — q;(£))v1 + 4:6'' 

j/3(e..o-i/3(e...)ii.< ^^^;:;^;/^^^., (3.14) 



+ — ^^^1=^ — '-^^<i 



-2, ,2 



VT+I? \/3(£:) -/i ^-a{e)J l + fi-e'^fi 

if e^fj, < 1. From (3.9), (3.13) and (3.14), for e > small enough and /i G (0, Aq) we 
obtain 

\\Hiw,e)\\^ < CiE)\\w\\^ + D{6), w G C^,o n 5,(0) (3.15) 

and 

\\H{w,€)-H{v,€)\\^<C{e)\\w-v\\^, w,v E C^^n B^{0), (3.16) 

where C{e), D{e) do not depend on /i and are given by 

C{e) = C^{e){l + M)+6+ \^V^^ - ^(") = ^'H^^'lUo- 

{/3{e) - Ao)Vl + 4£^ 

Since Ci{e) — )■ 0, f3{e) — > oo as £ — )■ 0+, by replacing 6 by 6/2 in (3.15), (3.16), we obtain 
(3.7) for e > sufficiently small. ■ 
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We now return to Eq. (3.5). Let < fi < Xq. For e > small, we look for a 
solution w G C^^ of (3.5). For the case yU = 0, where the space Cq^ denotes Cq, this 
question was addressed in [9]. Our purpose is to solve this problem for // G (0, Aq). 

We first apply a Liapunov-Schmidt reduction. From Lemmas 3.3 and 3.4, X^ : = 
Ker {J\c o) is finite dimensional, hence there is a complementary subspace F^ in C^,0; 

Cfj,fi = x^ © y^. 

For w G Cfj,fl, write w = C, + (p with C^ G X^,0 G Y^. Define S*^ := J|y^. Since 
S*^ : F^ — )■ C^fi is bounded and bijective, then S~^ is bounded. In the space C^,0) Eq 
(3.5) is equivalent to = S~^H{e,^ + 0), therefore we look for fixed points G F^ of 
the map 

J',{e,^,<P) = S;'HiE,^ + (/)). (3.17) 

For simplicity, in what follows we write S, J-", -Bcr(O) instead of S^, J-"^, -B^(O), respectively, 
when there is no risk of misunderstanding. 

Remark 3.7 For < /ii < /i2 < Aq with /ii,/i2 ^ 9f^o'(A), where o'(A) is the set of 
solutions of (2.2), it is clear that C^^ C C^^ with \\y\\^^ < ||l/||^2> ^'^'^ ^^^2 "^ ^/^i- 
Together with Lemmas 3.3 and 3.4, this implies that for each interval I := [fii, fJ^2] C 
(0, Ao) \ ^cr{A), we have X^^ = X^-^. We now show that the complementary suhspaces 
Y^ can he chosen so that Y^^ dY^ C. Y^-^ for ^ E I . In fact, let X^ = span {yi, . . . , yr}, 
where yi, . . . ,yr G C^^o (^''^'d r = r^ for fi E I . From the Hahn-Banach theorem, let 
hi G (C^i,o)' be such that hi{yi) = l,hi{yj) = for j ^ i,i,j = l,...,r. Define 
the natural injections i{fi,fii) : C^q ~^ C'^i,0; which are continuous, and the suhspaces 
Y^ = {y E C^fi : hi o i(/i,/ii)(|/) = 0,i = 1, . . . ,r}. Hence Y^ is a closed subspace of 
C^fl, and for y G C^,o we have XlLi hi{y)yi G X^,y - XlLi hi{y)yi G Y^, from which 
the decompositions C^^o = ^^i © Y^ follow, with Y^^ C F^ C F^^ for jj, E I . 

Theorem 3.8 Assume (H1)-(H4), and denote by (j{A) the set of characteristic values 
for (2.1). Fix an interval I := [fiifi2] C (0, Ao)\5Rcr(A), and denote r = r ^ for all fi E I . 
Then, there exist e* > and o" > 0, such that for < e < e* , the following holds: for 
each unit vector w eMP and all fi E I, in a neighbourhood Bj^{0) of u*{t) in C^, the set 
of all travelling wave solutions u{t,x) = il){ct + w ■ x) of (1.1) with speed c = 1/e and 
connecting to K forms a r-dimensional manifold (which does not depend on fi), with 
the profile ip E A^/,e, where 

Mi,e = {i^ : m =u*{t)+^ + 0(£, 0, for ^EX^n 5^(0)}, 

where (j){6, ^) = 0^(6, ^) is the fixed point of J^^{e, ^, •) in Y^ fl B^{0), and is continuous 
on (5,0- 

Proof In the sequel, we shall use the simplified notation S,J^,B„{0), for 
S^,J-'fj,,B^{0), respectively. Fix fi E I and k E (0,1). From Lemma 3.6 (cf. (3.15) 
and (3.16)), for S > small there are a = a{6,fi) > and e* = e*{6) > such that for 
0<e<e*,^EX^n 'Bjjyj and 0i, 02 e >^m ^ 5^(0) we have 

\\S-'H{^ + <j>,,e)\\ < \\S-'\\{C{e)U + M, + D{e))<S\\S-'\\{U + M, + l) (3-18) 
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and 

\\S-\H{^ + <P,,e)-H{^ + <P2,em,< ll^"lC'(e)||0i - 02lU < ^ll-^-l ||0i - 02lU(3.19) 

with 5(1 + 2a)\\S-m < a and S\\S-'^\\ < k. From (3.18) and (3.19), it follows that 
J" : (0,£*) X (X^ n -Bct(O)) X {Y^ n 5^(0)) -^ Y^n B^{0) is a uniform contraction map 
of G F^ n B^{0), hence for (e,^ ^ (0)^*) ^ (-^^ l~l B^{0)) there is a unique solution 
(p{e,^) = (f)^{e,^) G Yfj_ of (3.17), with 0(6,0 continuous. Define the r-dimensional 
manifold M^^s = {^ ■ ^{t) = u*{t) + ^ + 0(5,0, for ^ G X^ n 5^(0)}. Choose 
0" = a{6,fi2), independent of yU G /. From the uniqueness of the fixed point and 
Remark 3.7, it follows that 0^(^,0 = 0/^2(^)0 does not depend on /i G /, as well 
as M^^s ■= Mi^e- ■ 

We observe that if is a hyperbolic equilibrium of (1.2) and / has the particular 
form /(0) = F{(j){0),g{N(j))), then Theorem 3.1 asserts that the result in Theorem 3.8 
is valid for yU = 0. 

Corollary 3.9 Under the assumptions of Theorem 3.8 and with the same notation, for 
< /i < Aq such that the strip {A G C : 3?A G (;U, Aq)} does not intersect (y{A), the 
manifold Aifj,^e is 1- dimensional. 

Corollary 3.10 Under the assumptions of Theorem 3.8 and with the same notation, 
for an interval I := [/ii/U2] C (0, Aq) \ 3?(t(A), there are e* > 0, o" > and C > such 
that the travelling profiles ip{e,C,) satisfy 



{e,Oh<C, \\^'{e,Oh<C for < e < e*,^ e X,n B^{0), 
where C does not depend on fi E I. In particular \ip{e, C,){t)\ < Ce^* fort < 0,0 < e < e*, 

^eX^nBjo). 

Proof For all fi & I, the profiles are given by -0(5,0 = u* + ^ + 0(5,0, where 
0(5,0 = 0/.2(^>O is the fixed point of ^'^2(5,^, ■). Since \\y\\^ < \\y\\^^ for y G 0^,^^, we 
only need to prove the result for fi = H2- In what follows, we write ST^^ = S^^. 

Fix k G (0, 1), and consider Si > such that ||S'^"'^||C(5) < k for < 5 < 5i where 
C{e) is as in (3.16). From (3.19), for wi,W2 G C^2,o H -Bo-(O) we have 



IS" ^{H{wi,e) - H{w2,e))\\^2 < k\\wi - W2\ 



M2! 



and the contraction principle yields 

i^ll^..(^,e,o)IU2<Y^^ 

For 5, 0" > small enough, from (3.15) we get 



e,Oh. < T^Tll'^^2(^,e, 0)11^2 < t^tII^"'IIII^(^'OIU., e e X,nB^{0). 



\\Hie,Om,.<Cie)m,,+e'\\u*"\\^„ ^eX^nB,{0). 

We thus obtain 11-0(5,011^2 ^ ll'"*IU2 + i'^ + C {e))cr + e'^\\u*" \\^^ < Ci for 5 small, where 
Ci does not depend on 5, C,- 
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Now we want to prove a similar estimate for the derivates dilj{e,^)/dt. For 
simplicity, we only prove the result for .^ = 0. 

Since tp(t) := ■ip{e,0)(t) is a solution of (3.2), then ip(t) is given by the integral 
formula 



from which we derive 

t +00 



^'(t) = A^ {a{e)l e"(^)(*-)[V'(s) + /(^.)]rf. - ml e^^^^^^^-^H^l^) + f{i^s)]ds). 



— oo i 

Since / is bounded on bounded sets of C = C{\—t^ 0]; M^), there is £ such that |/(0)| < d- 
for e C with ||0||oo < C*!- Thus, \f{;ips)\ < £ for s G M, where £ does not depend on yu, e, 
and IIV-'lloo < 2(Ci + £)/Vl + 4e2. From (3.10), the C^-smoothness of / and /(O) = 0, 
we easily deduce that there is C2 > such that \\ip'\\^ < C^- This completes the proof." 
In fact a stronger result can be proven: 

Corollary 3.11 Assume (H1)-(H4), take fi E I := [fJ,ifJ,2] C (0, Aq) \^a{A), and 
consider the travelling wave profiles ip{e,^) = u* + ^ + (f){e,^) for e G (0,e*),^ G 
X^ n -Bcr(O), given in Theorem 3.8. For .^ = and ^ E I , the profile ip{e, 0) satisfies 

ipie, 0) -^ u* in C^ as e ^ 0+. (3.20) 

Proof. Let e*,^ > be as in the statement of Theorem 3.8, and recall that 
ip{e, 0) = ipn{e, 0) only depends on e. Next, we deduce some estimates as in Lemma 3.6, 
so details are omitted. For e = 0, define 

H{0,w){t)= [ e-^'-^\f{ws + K) - fiK) - Df{u:)ws]ds. 



We write H{0,w){t) = f^^e-^^-'^G{0,t,w) ds := H2{0,w), where G{0,t,w) is given 
by (3.4). After some computations, we observe that the function H restricted to 
[0,e*) X {C^fiDB^iO)) satisfies 

\\H{e,w) - H{0,w)\\, < Co{e)\\w\\, + Do{e), 

with Co{e),DQ{e) independent of /i, Co{e),DQ{e) — )> as £ — )> 0+. This means that the 
function {e,w) t-)- H{e,w) converges, uniformly on tf; G C^,o H -B^(O), to H{0, ■) in C^^o 
as £ — )■ 0+. 

Moreover, for e = and ^ = the fixed point of (3.17) is 0(0, 0) = 0. Therefore, the 
application of the contraction principle as in the proof of Theorem 3.8 leads to (3.20)." 

Remark 3.12 As seen in Section 2, the existence of a positive eigenvector v G M^ 
associated with the characteristic root Aq of (2.2) was crucial to prove the existence of 
a positive heteroclinic solution u*{t) of (1.2), connecting the equilibria to K . For all 
the results in this section the positiveness of such heteroclinic solution is irrelevant, and 
therefore it is not necessary to impose the above requirement in (H4) that v is positive. 
For the same reason, in Section 3 assumption (H2)(ii) is not needed as well. 
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4. Positiveness of travelling waves 

Consider tlie cliaracteristic equation for the linearization of (3.2) at 0, 

det A,{z) = 0, where A,{z) := e^z^I - zl + L{e^'I), (4.1) 

where L = Df{0). Recall that for e = the characteristic matrix- valued function 
Ao[z) was defined in (2.2). Since Aq > is a simple root of the characteristic equation 
det Ao{z) = 0, from the implicit function theorem, for e > small there is a simple real 
root \{e) of (4.1), with \{e) -)■ Aq as £ -)■ 0+. 

Lemma 4.1 For 6 > sufficiently small and Si > 0, there exists ^o > such that, for 
< e < Eo, X{s) is the only root of the characteristic equation (4-1) on the vertical strip 
\q — 5 < Re z < \q + 5i. 

Proof. Let 5 > be such that Aq is the only root of det Ao(-2) = on the strip 
S = {z : \q — 5 <^z < \q + 5i}. If z{e) G S* is a root of det A^{z) = 0, then there is a 
unit vector w = w(e) G M^ such that {e'^z^ey — z(€))w = L(e^^'^^'w), hence 

\\L\\ > \e^z{ef-z{e)\ > \'^ {eh{ey - z{e))\ = \'^ z{e)\\2e^^z{e) - 1\. 

Choose £0 > such that |2e^5R2; — 1| > 1/2 for all z E S, \e\ < Eq. For |e| < Eq, we have 

\'^z{e)\ <2\\L\\. 

Thus, for \e\ < Eq the solutions z{e) G S* of detAi;{z) = are necessarily inside the 
rectangle F = [Xq ~ 6, Xq + 5i] x [ - 2||L||, 2||L||]. 

Now, let F{z,e) := detAe{z),z G C,e G M. Clearly, F{z,e) -^ F{z,0) as £ ^ 0, 
for all z G C Moreover, since Ai,{z) = e'^z^I + Ao(-2), one easily deduces that the 
function F(-, e) converges uniformly to F(-, 0) on bounded sets of C, as e — )■ 0. 

We now apply Rouche's Theorem on the boundary dV of F. Set m = 
min^gar |-^(-2, 0)1 > 0. For \e\ small, we have \F{z,e) — F{z,0)\ < m, z E dT, hence 
F{z, e) and F{z, 0) have the same number of zeros inside F. Thus, for e > sufficiently 
small X{e) is the only solution of (4.1) in the strip S. ■ 

For (3.2) written as a system in IR , its linearized equation at zero is 

x'{t) = Ce{Xt), (4.2) 

where 



For the linear system (4.2), the characteristic equation is given by 

det aw = 0. where DM={^^\i,.,^ (s -4,) ^ ' '"' 

and I is the N x N identity matrix. Clearly, det Ai;{s) = e'^^ det D^{s), hence for e > 
(4.3) is equivalent to (4.1). 
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Lemma 4.2 Consider b E M. and Eq E (0, 1) such that det Aj(s) ^ for all e E [0,£o] 
and s on the vertical line H = {s = b + iy : y E M}. Then there is e\ E (0, Sq] such that 

sup < \s\ ||A£(s)"^|| : < e < ei, s G s| < oo, 



and 



sup I |s| ll-De(s) ^11 : < e < £i,s G sl < cx). 



Proof. For e G (0, eo],s E S, we have 



«(-)---^«-("r ..voc-V/';^::.)' 



^o(s)-' \ / -/ 0\ __^,^ , „. _n+ 



Clearly, 

''<"'- r"o' AoW-JUiCe-/) oJ^=''<«-' ^" ^^« 

with G(e, s) continuous on [0,£o] ^ S- For s G S, we have ||L(e'^'/)|| < max(l, e~^^)||L||. 

It follows that 

\\Gi6,s)\\ < cie%s\ + l)||A,(s)~l, < £ < 5o,s e S, 

for some c > 0. Since Ae{s) = (e^s^ — s)I + L{e^'I), then 

"^'^^^ " " |s|k2s-l|-||L(e-J)|| 

if Islle^s - 1| - ||L(e^7)|| > 0. Choose si E (0,£o] such that 1 - ejb > 1/2. 
Then for s G S and < e < ei, if \s\ > 4max(l, e^^^)||L|| := ci, it follows that 
\s\\e^s - 1| - ||^(e"-/)|| > |s|/4 > 0, thus \s\ \\A,{s)-^\\ < 4 and 

I I 11^/ \ii ^ cIsKe^lsl + 1) ^ ce^IsP 



for e E (0,ei] and s eTj, \s\ > Ci. Now, \s\\e^s - 1| > |s|a/£:''|sP + 1/2 > e^jsp, and 

<2 
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|s||£2s-l|-||L(e^-J)|| 
if £2|s|2 > 2||L(e^7)||; andif e^jsp < 2||L(e^7)||, then 



|s||£2s-l|- ||L(e^7)|| - |s|/4 

hence, |s| ||G(e, s)|| < 6c if |s| > Ci. 

On the other hand, on the compact set {{e,s) E [0,ei] x S : |s| < ci} the 
continuous functions |s| ||A£(s)^-'^|| and |s| [[^(e, s)|| attain their suprema, and the 
conclusion follows. ■ 

We are finally in a position to prove the main result of this section, on the existence 
of positive travelling wave solutions of Eq. (1.1) for large wave speeds. 
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Theorem 4.3 Assume (H1)-(H4)- Then, there is c* > 0, such that for c > c* Eq. (1.1) 
has a positive travelling wave solution of the form u{t,x) = il){ct + w ■ x) for each unit 
vector w G W, with ■?/'(— oo) = O,ip{oo) = K. Moreover, the components of the profile 
ijj are increasing in the vicinity of —oo and it satisfies ^{t) = 0(e'^*^^^*),-?/''(t) = 0(e'^*^^^*) 
at —oo, where e = 1/c and \{e) is the real solution of (4-1) with \{e) — ?■ Aq as e — )■ 0^. 

Proof. Consider Eq. (3.2), where e = 1/c. Let fi G (0, Aq) be as in the statement of 
Corollary 3.9 and satisfy Aq < fi + 26 for some fixed 6 G (0,/i/4). Suppose also that X{e) 
is the unique solution of det Air{z) = on the strip fi — 6 < ^ z < 2fi for all e E [0, e*]. 
Fix the profiles 

^,=^(£,0), ^o(t):=M*(t) = e"«*v + 0(e2^*), ee[0,e*], (4.5) 

as in Corollary 3.11 and Theorem 2.1. Recall that v is a positive eigenvector associated 
with Aq. The proof is now divided in several steps. 

Claim 1. There is eo > such that (^'^(t), V^^(t)) = (e^(^)*vi(£), A(e)e^(^)*vi(£)) + 
We(t), e G [0,eo), with continuous \'^{e) > 0, v^(0) = v, and We{t) = 0(e(^"+'^)*) at -cx). 

To prove the above claim, note that x^it) := (■^/'^(t), ■?/'^(t)) is a solution of the system 

x[(t) = X2(t) 

e'x'^it) = X2{t) - L{xi^t) - heit), (4.6) 

where L = Df{0) and he(t) = f{{ipe)t) — L{{ips)t)- For e > 0, equivalently we write 
(4.6) as 

1 / 



x'{t)-C,{x,) ^2 ^^^^^^ 

where C^ is as in (4.2). Since / is a C^ function, using the Taylor formula for / (cf. e.g. 
[4, p. 23]), we have the estimate 

\he{t)\< f il-s)\\D'f{sii,,)t)\\\\{4,,)t\\lds, teR,ee[0,e*]. 
Jo 

Since Wipe — M*||/i — )■ as e — )■ O"*" and u*(t) — )■ as t — )■ — oo, the continuity of V^f at 
implies that ||D^/(s(-?/'e)t)|| is uniformly bounded on 5 G [0,ei] C [0,£*), t < 0, s e [0,1]. 
Together with Corollary 3.10, this leads to 

\xe{t)\ < Ce^"', \K{t)\<De^^' for t < 0, (4.7) 

for some constants C, D independent of e G [0, ei\. 

We now apply Theorem 6.2 to (4.6) at — oo (see Appendix), and derive that for 
£> 

Xe{t) = {iJe{t),ij',{t)) = ^e{t) + We{t) (4.8) 

where Zs{t) is an eigenfunction for the linear system x'{t) = C^{xt) corresponding to the 

set Ke = {z eC: det A£(^) = 0,/i < 3?^^ < 2;u} and Weit) = 0(e(^"+'^))*) at -oo. From 



Positive travelling fronts for delayed systems 18 

Lemma 4.1, let e be on some interval (0,£:o) C (0,ei) such that that A^ = {A(e)}. Then, 
Zeit) is an eigenfunction for x'{t) = Ce{xt) associated with the root \{e) of (4.3), hence 
z,{t) = e^(")*v(e) with v(e) = {v\e),Y\e)) G M^a^ satisfying D,(A(£))v(£) = for D, 
as in (4.3). From this we obtain v^(e) = X{e)\^{6) and A^{X{e))\^{e) = 0. Furthermore, 
from Theorem 6.2 (with a = fi,b = \o + 26,e = 6) and formulae (6.5) and (6.7) (adapted 
to the situation — oo) we get 

Ze{t) = Res (e* x„ X{e)) = e^^'^' lim (s - X{e))Si,{s) = e^^'^^Y^e), X{e)Y\e)) 

s— !>A(e) 



and 



with 



b+4+Joo 



Mt) = 7^ [ ' e'%is)ds, (4.9) 



Ms) ■■= Xsi-K-s) = G{e, s) { r,{s) - ( ~ ° ) ) , (4.10) 



'■e\ 



for G{€,s) = Di;(s) ^ and 



r,{s) = xM + A e'- / e-^"x,(M) rfw , 











hs{s) := hs{—){-s) = e~'''he{u)du, ^s<Xo + 26. (4.11) 

J —oo 

Here, X£( — ),hs{ — ) denote the Laplace transforms of the functions t h-)- Xi,{—t),t i— )■ 
hi;{~t), respectively. Note that Xe(s) is meromorphic for ^s < 2fi with a unique 
singularity at s = X{e) which is a simple pole of As;{s)^^ (cf. Appendix). 
From (4.8), we get 

Mt) = e'^^K\e) + wlit), ij'M = X{e)e'^^K\e) + w',{t), (4.12) 

with wl{t) = 0(e(^«+'^))*) at -oo and w^{t) = (wKt))'. Next, the definition of r^{s) in 
(4.11) yields 

rei-s) =3:^(0)+ ( 1 ./ „. .0 



\ 



hence from (4.4) and (4.10) we obtain 

/ A,{sy'L(e'fe-''Mt)dt] \ / ^ fs^ih (s) \ 

Us) = G{e,s)xM - V n {]-{ : \ iZ \ ) ■ (4-13) 

\sA,{s)-^L(e-fe-^'Mt)dt)j \sA,is)-^h,is) J 

We now extend naturally this situation for e = 0. In (4.5), denote A(0) = Aq, v^(0) = v. 
Write Xe = (x^,x^) for e G (0,eo); and let ^l{s) be defined by (4.13) for e = 0. Note 
that formula (4.10) can still be used to obtain Si.l{s) (cf. Appendix for more details), 

3cl{s) = Aois)-'His) + Ms)], 

where Vq^s) = u*{0) + L^e"' f e~''^u*(t) dt) = lim£^o+ rl{s). 
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For each e G [0,£:o)) s = \{e) is a pole of order one of G{e,s), and from (4.10) 
we deduce that for e G [0,£o) and /i < $Rs < Aq + 25 the function A{e,s) defined by 
A{e, s) = (s — A(e))x^(s) for s ^ A(£:), A{e, \{e)) = v^(e) is analytic on s and continuous 
on (e, s). In particular, \im.!r^Q+ A{e , \{e)) = A{0,\q) = v > 0, hence ^r^{e) — )> v as 
e — )■ 0+. Moreover, v^(e) > for e > sufficiently small. This proves Claim 1. 

Claim 2. For ^o > sufficiently small, there exists a constant Dq > such that 

\wl{t)\ < Doe(^»+'^)* for all t < 0, < e < £*. (4.14) 

To prove Claim 2, once more we shall use some formulae and estimates in the proof 
of Theorem 6.2 in the Appendix, changed accordingly to account for the asymptotic 
behaviour at — oo, rather than oo. 

Define Ve{t) = {vl{t),v^^{t)) = e-(^«+'^)*«;,(t) and M,(t) = {ul{t),u'i{t)) = 
g-(Ao+35/2)t^^(^)^ Note that vl{0) = wl{0) an vl{t) = wl{0) - J^{viy{s)ds, with 
wl{0) = ^-£(0) - ^r\e). Then 1^^0)1 < C + \^r\e)\, where C > is as in (4.7). 
We need to prove that vl{t) is uniformly bounded for t < and e > small enough. 
In order to achieve this, we shall show that there are constants Co, Dq > and ^q > 0, 
such that 

lk'l|Li(-oo,o] <C'o, ee{0,el) (4.15) 

and 

||(^^.')lUi(-oo,o] < Do/2, e G (0,6*), (4.16) 

so that (4.14) follows immediately from (4.16) and |w^(0)| < Dq/2 for e G (0,£o). These 
uniform estimates require a careful analysis of the explicit formulae for We given in (4.9) 
and (4.10). We shall prove (4.15) beforehand, and then use (4.15) to prove (4.16). 

First, observe that a;e(t) = -2e(t) + '^^(t) is a solution of (4.6), with Zi.{t) being an 
eigenfunction for the linear system x'{t) = Ce{xt), hence We{t) is a solution of system 
(4.6) as well. The definition of v^{t) yields now 

{vl)\t) = -{\o + 6)vl{t)+vl{t) (4.17) 

and 

e\vlnt)-a{vly{t) + PM = ^, (4.18) 

where a = 1 — 2£^(Ao + 5) and 

P,(t) = [e\\o + 5f - (Ao + 5)]vl{t) + L(e(^«+^)-(^')*) + e-'^^^+'^'Kit)- (4.19) 

Next, similarly to (6.9) and Remark 6.3, using (4.9) and the Plancherel theorem, 
we obtain 

Jt/2 r+oo 



w 



£ I 



^i^^^^_(Ao+5)<||^ = 11^— / e^"*5,^(Ao + 35/2 + lu) du\\ 



Li(-oo,0] = WeV')^ llLi(-oo,0] = 11^— / « X^ I^Aq "h 00/ ^ "h 'mj a'U||Li(_oo,0] 

<Ci(5)||X^(Ao + 35/2-2-)|U2(iR) for £G(0,£o). (4.20) 
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From (4.7), |a;e(0)| is uniformly bounded for e G (0, Eq); and since 2/i > \q + 2S, (4.7) also 
implies that ||L(e^' J e~^Ve(^) '^M II ^iid |h£(s)| are uniformly bounded for e E (0,eo) 
and 3?s = Ao + 35/2. Using now Lemma 4.2, from (4.13) we conclude that there is 
Ki > such that 

or 

|x^(s)| < Ki/\s\, ior s = b + — +iy, y G M, and < £ < eo, 
from which we derive 

P^(Ao + 3V2-^-)lli^(iR)< / f.^^L^ , dy<oo. 

Jr (o + fr + y^ 

Together with (4.20), this leads to the estimate (4.15) with ^q = Eq. 

In order to prove (4.16), we now use (4.18) and some ideas from Aguerrea et al. [1, 
Lemma 4.1]. After partial integration of Eq. (4.18), we obtain 

{viy{t) = e(°/^')*(t;,i)'(0) + —,— / e-("/^')^P,(s) ds. (4.21) 



^ Jt 



Let e > be small, so that a > 0. From (4.17) and Claim 1, 

{vl)'{0) = -(Ao + 5)wl{0) + wliO) = -(Ao + 5)^0) + {\o + 5- \ie))^^\e) + ^^(0), 

with v^(e) — )■ V as e — )■ 0+. Corollary 3.10 implies that |('i^£)'(0)| is uniformly bounded 
on e > sufficiently small, hence 

liviym f e("/^')*rft<ir2, (4.22) 

for some K2 > and e > sufficiently small. On the other hand, interchanging the 
order of integration leads to 

J_ / gW^')* f f e-^''/''>Pe{s)ds\ dt = \ j e-(°/^')^P,(s) ( T e^"/^')* rft") ds, 

hence 

J_ f (,(^/e^)t ( I g-(a/.2)s|p^(^)| ^^ dt<- I \Pe{s)\ds. 



— 00 \J t 



From (4.7), (4.15) and the definition of Pe{s) in (4.19) , we easily see that 

\\Pe\\LH-oofi] < K3, <e <el 

for some ^q > 0. Together with (4.21) and (4.22), this yields the estimate (4.16), and 
therefore Claim 2 is proven. 

We finally prove: 

Claim 3. There is e* > such that ipeit) > for t G M and e G (0, £*). 

From Claims 1 and 2, for £ > small enough 
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where 1 = (1, . . . , 1) and |A(e) - Ao| < S/2. Choose T* < and e* > such that 
Y\e) > v/2 if < £ < el, and Doe^^^^^^"! < v/4. 

Then 

Ait) > e^^'K/4 > 0, 0<e <el, t<T*. 

On the other hand, since Ij-i/'e — ■u*||oo — ^ as £ — )> 0"^, we define rj := mit>T* u*(t) > 0, 
and suppose that €{ was chosen so that \\ipe — u*\\oo < 1] ior < e < el- It follows that 
ipeit) > for all t G M and e G (0, e*). The proof of the theorem is complete. ■ 

The above proof shows that the requirement that there is a positive eigenvector 
V G M^ for the dominant characteristic value Aq of (2.1) is crucial to deduce the 
positiveness of the travelling wave fronts, for large wave speeds. Nevertheless, the 
existence of such waves and their asymptotic behaviour at — cxd can be deduced from 
our Theorem 3.8, as well as the auxiliary results in Section 3, and the proof of the above 
Claim 1. We summarize these remarks in the following theorem: 

Theorem 4.4 Assume (HI), (H2)(i), and 

(i) for Eq. (1.2) the equilibrium u = K is localy asymptotically stable; 

(a) for Eq. (1-2), the linearized equation about has a real characteristic root 
Aq > 0, which is simple and dominant; 

(Hi) Eq. (1.2) has a heteroclinic solution u*{t),t G M, with u*{—oo) = 0,u{oo) = K 
andu*{t) = 0(e^o*) at -oo. 

Then, there is c* > 0, such that for c > c* , Eq. (1.1) has a travelling wave 
solution of the form u{t,x) = ilj{ct + w ■ x) for each unit vector w G W\ with 
^{t) = 0(e'^*^^^*),-(/''(t) = 0{e'^^^^*) at — oo, where e = 1/c and A(e) is the real solution of 
(4-1) with \{e) -> Ao as e ^ 0+. 

5. Applications 

5.1. A diffusive generalized logistic equation with distributed delay 

As a first application, we consider a scalar reaction-diffusion equation with distributed 
delays in the reaction-terms, which includes the Fisher-KPP equation with delay as a 
particular case. 

Let C = C{[~T, 0]; M), r > 0, and consider 

Ou Ou 

-^{t,x) = —{t,x)+bu{t,x)[l-L{ut{-,x))], teR,xeR, (5.1) 

where 6 > and L : C — > M is a nonzero positive linear operator, i.e., L 7^ is linear 
and L{ip) > whenever v^ > 0. In particular, L is bounded and ||L|| = L(l) > 0. 
The corresponding delayed ODE model, 

u{t) = bu{t)[l - Lut], t G M, (5.2) 

has two equilibria, u = and u = L(l)"^ := K. Here / in (1.2) reads as 
/((/?) = bip{0)[l — L{ip)],ip G C, for which it is easy to verify that conditions (HI) 
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and (H2) are satisfied. The linearized equation about zero is the ODE u'{t) = bu{t), 
with characteristic value 6 > 0. 

Now, u'{t) = —bKLut is the linearized equation about K, with characteristic 
equation P(A) := A + hKL{e^') = 0. If A = ioj.oj > 0, is a solution of -P(A) = 0, 
then L(cos(w-)) = and = w + hKL{sni{uj-)) > u — b. If ur < 7r/2, we deduce that 
there is 5 > such that cos{uj9) > 6 for 9 E [— r, 0], hence L{cos{uj-)) > SL{1) > 0, 
which is a contradiction. On the other hand, if br < 3/2, from [8] we conclude that the 
positive equilibrium K is globally attractive in the set of all positive solutions of (5.2) 
with initial conditions (/? G C_|_, V'(O) > 0. We thus conclude that (H3) holds if br < 3/2. 
Therefore, the following result is an immediate consequence of Theorem 4.3 

Theorem 5.1 If br < 3/2, there exists c* > such that for c > c* equation (5.1) has 
a positive travelling wave solution u{t,x) = iIj{x + ct) with ■?/'(— oo) = 0,-?/'(oo) = K. 
Moreover, ip is increasing in a vicinity of —oo and it has the asymptotic decay il){t) = 
0(e''(^)*),V^'(t) = 0(e^(^)*) at -cx), where 6(c) = 26c/(c + Vc^ - 46). 

We note that (5.1) includes as a particular case the Fisher-KPP equation with a 
single delay, 

Ou Ou 

{t,x) = —^(t,x)+bu(t,x)\l-u(t-T,x)/K], teR,xeR. (5.3) 
ot ox^ 

By using a pair of upper-lower solutions and a monotone iterative method, Wu and Zou 

[23] proved that if c > 2\/b, then there exists t*{c) > such that, for a delay r < t*{c), 

(5.3) has a non-decreasing travelling wave front connecting to -ft' with wave speed 

c. Our approach does not allow us to determine the minimal wave speed c*, but, on 

the contrary, we explicitly exhibit the maximal delay r* = 3/(26), under which we can 

assure the existence of such positive (but not necessarily monotone) travelling solutions. 

Travelling waves for the Fisher-KPP equation (5.3) with 6=1 were also considered 
in [9, Corollary 6.6], where it was shown that, for r < e^^, there exists c* > such that 
there is a travelling wave front with wave speed c > c*. We remark that Theorem 5.1 
applied to (5.3) with 6=1 clearly improves this result, since it guarantees the existence 
of travelling wave solutions for r < 3/2, and, most relevant in biological terms, it does 
assert that such travelling waves are positive. 

For some other recent results and references on the Fisher-KPP equation, see [2, 12]. 

5.2. A chemostat model with delayed growth response 

Consider the following model for the growth of bacteria in a well-stirred chemostat 
supplied by a single essential nutrient (cf. Ellermeyer [6] and Ellermeyer et al. [7]): 

S'{t)=D{S'-S{t))-f{S{t))u{t), 

u'it) = e-^^f{S{t - T))u{t - r) - Du{t). (5.4) 

Here S{t) and u{t) are the concentration of nutrient in the growth vessel and the biomass 
concentration of bacteria at time t, respectively, -D > is the dilution rate of the 
chemostat, 5*° > is the input concentration of nutrient, and r > is the delay in the 
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growth response, to account for the lag in the nutrient conversion into biomass due to 
cellular absorption; / is the specific functional response for the bacteria, and typically 
the Michaelis-Menten response is chosen, f{s) = ms/{a + s), s > 0, with m,a > 0. 
More generally, one can consider a continuously differentiable and bounded function 

/ : [0, oo) -^ [0, oo) with 

/(O) = 0, f{s) > for s G M. (5.5) 

For an unstirred chemostat, the nutrient is added to the vessel but not mixed, so one 
has to introduce diffusion terms. The diffusion rates di,d2 > for the nutrient and the 
organisms may be different, and model (5.4) becomes 

— (t, x) = diAS{t, x) + D{S^ - Sit, x)) - f{S{t, x))u{t, x) 

-_(t, x) = d2Au(t, x) + e-^^fiSit - T, x))u(t -t,x)- Duit, x), (5.6) 
at 

for t G M and x G (0, L) {L > 0) (or more generally x & Q, where fi C M^ is an open 

domain). There is an extensive literature on chemostat models with "delayed growth 

response", with and without diffusion. We refer to [3, 6, 7, 11, 21, 22], for results, other 

related chemostat models, biological explanations, and further references. 

For both (5.4) and (5.6), there is always the equilibrium (S"^,0), corresponding to 

the "washout state". If 

f{S') > De^\ (5.7) 

there is another nonnegative equilibrium (S*, -u), called the "survival state", given by 
{S,u) = {f^^{De^'^),e~^'^{S^ — S)). Condition (5.7) imposes a restriction on the size 
of the time-delay r, which should satisfy r < D^^ \og{ f{S^)D^^) for / as in (5.5). 
Moreover, (5.7) implies that the equilibrium (5*°, 0) of (5.4) is unstable and {S,u) 
is asymptotically stable and a global attractor of all solutions with initial conditions 
(5o,Wo) = (0i,02)gC+, 02(O)>O[6, 7]. 

In order to apply our results, we first observe that both the positive cone C+ and the 
set {(01,02) £ C+ : 02(0) > O,0i(O) < S'^} are positively invariant for (5.4). Translating 
the washout state to the origin, by setting s(t) = S^ — S(t), we rewrite (5.4) as 

s'{t) = -Ds{t) + /(5° - s{t))u{t) 

u'it) = e-^^f{S° - s{t))u{t - r) - Du{t). (5.8) 

If / is C^-smooth and (5.7) holds, then Eq. (5.8) satisfies (HI) and (H2). The set 
A = {(01,02) G C+ : 02(0) > 0,0 < 0i(O) < S^} is positively invariant, and (5.8) 
has equilibria Eq = (0,0) and K := {s,u) > 0, with s = S^ — S,u = e~^^s, the first 
one being unstable and the second one being locally stable and a global attractor of all 
solutions with initial conditions in A (cf. [6]). It remains to verify that (114) holds. 
The characteristic equation for the linearization of (5.8) at (0, 0) is 

(A + D){X + D- e-^V('5°)e^^") = 0. (5.9) 

Define h{x) = (x + D)e^^'^^^^,x G M. Under (5.7), there are two real roots of (5.9), —D 
and Ao, where Aq > is the unique solution of h{x) = /(S"^). For A ^ M a solution of 
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(5.9), we have h{^\) < f{S^). Since h'{x) > for a; > 0, if follows that 3fJA < Aq, and 
we conclude that Aq is a dominant eigenvalue for the linearization of (5.8) about (0,0), 
with V = (/(S'*')(Ao + D)"^, 1) > as an associated eigenvector. From Theorems 2.1 
and 4.3, we therefore obtain the following result: 

Theorem 5.2 Consider Eq. (5.4), where S^, D,t > 0, the function f : [0, oo) — )■ [0, oo) 
is bounded, C"^ -smooth and satisfies (5.5) and (5.7). Then, there exists a heteroclinic 
solution {S* (t) , u* (t)) connecting the washout state (S'*',0) to the survival state {S,u) of 
(5.4), with < S*{t) < S^,u*{t) >0 forteR and {S* (t) , u* (t)) = (^°,0) + 0(e^o*) at 
— oo, where Aq > satisfies \o + D = /(S'°)e^*^'^°+^''^. For the diffusion model (5.6) with 
di,d2 > 0, there is c* > such that for c > c* (5.8) has a positive travelling wave solution 
of the form {S {t , x) , u(t , x)) = {ipi{ct + x) , ■i/j2{ct + x)) , with ■?/'i(— oo) = S^,ip2{~oo) = 
and ipi{oo) = S,ip2{oo) = u; moreover, ip'iit) < Q^ip^it) > in the vicinity of —oo, and 
{ipi{t),ip2i't)) = ('S'*',0) + 0(e^*^'^''*) at — oo, where z(c) is the real solution of 

z^-cz- c{D - e-(^+^)V('5°)) = 0. (5.10) 

Proof. Let e = 1/c, for c > large. With the notation in (4.3), we have that 

det DJX) = — l-(5rfiA' -X-D) (eX^ -X~D + e^(^+^)V('5°)') . (5.11) 

Define X{e) as the real solution of (5.11) such that X{e) — )■ Aq as £ — )> 0+. Then 
z{c) = X{e) satisfies (5.10). ■ 

Theorem 5.2 asserts the existence of positive travelling wavefronts for (5.6) with 
S{t,x) < S*" for all t G IR, a; G IR. Here, due to the change of variables s = S — s^, 
the positivity of the component s{t) (or s{t,x)) translates as the nutrient concentration 
being smaller than S^. We emphasize that biologically significant solutions of (5.4) and 
(5.6) must be positive and have a nutrient concentration S not larger than the input 
concentration S"". 

Remark 5.3 The existence of a positive eigenvector v associated with the dominant 
eigenvalue Xq, as prescribed in (H4), may seem a quite restrictive requirement, since it 
is not satisfied by many populations dynamics systems, namely Kolmogorov type models 
with N > 1. We however observe that if the characteristic matrix for (2.1) at Aq, 
Ao(Ao), is an irreducible matrix with non-negative off-diagonal entries, then there is a 
positive eigenvector for Aq{Xo) associated with Xq (see e.g. [21, p. 258]). This property 
will be exploited in a forthcoming paper, where Theorems 3.8 and 4-3 will be applied to 
several population models. 

6. Appendix 

In this appendix, we extend Proposition 7.1 of Mallet-Paret [18] to systems with 
distributed delays. 

Consider the FDE 

x{t) = LQXt + h{t), teR (6.1) 
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and the homogeneous system 

x'{t) = Loxt, (6.2) 

where Lq : C([— r, 0]; M^) — )■ M^ is a bounded hnear operator and /i : M — )■ M^ is 
continuous. For (6.2), write the characteristic equation 

det Ao(s) = 0, where Ao(s) = si - Lo(e'7). 

It is well known that the solutions of the characteristic equation are exactly the 
eigenvalues for the homogeneous system (6.2), i.e., the eigenvalues for the infinitesimal 
generator A associated with the semiflow of (6.2). Furthermore, the spectrum (t{A) of 
A is only composed of the point spectrum. 

Lemma 6.1 If f is a holomorphic function on a disc {s : \s — \\ < e}, where A is an 
eigenvalue of (6.2) and e > zs small, then 

x{t) = Res (e*-Ao V, A) = x^ / e^*Ao(s)- V(s) ds 

2vr« J\s-x\=e 

is an eigenf unction of (6.2) corresponding to A. 

Proof. The proof follows the arguments of Mallet-Paret [18, Section 7], so we omit it. ■ 

Theorem 6.2 Let x{t) be a solution of (6.1) on [T, oo) for some T eM.. Assume there 
are a,b EM.,a < b, such that 

x{t) = 0{e-'''), h{t) = 0{e-^') as t ^ oo. 

Then, for every e > 0, we have 

x{t) = z{t) + 0{e-^^^'^') as t ^ oo, 

where z{t) is an eigenfunction of (6.2) associated with the set of eigenvalues A = {A G 
cr{A) : —b <^\ < —a}. Analogously, if x{t) is a solution of (6.1) on (— oo,T] for some 
T eR and 

x{t) = 0(e"*), h{t) = 0{e''') as t -^ -oo, 

with a < b, then for every e > we have 

x{t) = z{t) + 0(e(^~')*) as t -^ -oo, 

where z{t) is an eigenfunction of (6.2) associated with the set of eigenvalues A = {A G 
a{A) :a< 3f^A<6}. 

Proof. We only prove the result for +oo; for — oo it is analogous. Without loss of 
generality, take T = 0. In what follows, for / : [0, oo) — )> C,f{t) = 0(e^"*) at +oo, we 
denote the Laplace transform of / by 

_ /"OO 

(£/)(s) = /(s)= / e^''f{t)dt, for ?fis>-a. 
Jo 
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Write Lq{ip) = J dri{9)ip{9), where r]{6) is an N x N matrix- valued function of 
bounded variation. Applying the Laplace transform to (6.1), we obtain 

-a;(0) + sx{s) = C{LoXt){s) + h{s), ^s> -a, (6.3) 

where 

C{LoXt)is) = / e-'' i / dri{9)x{t + 9)\dt= / dri{9) i / e-"(*+^)x(t + 9) dt 

d7]{9) ( f e-'''x{u)du + x{syj =Lo(e'- f e"''x{u)du] + Lo(e^7)£(s). (6.4) 
From (6.3), (6.4), we obtain 

Ao{s)x{s) = r{s) + h{s), where r{s) = x{0) + Lq [ e^' e'^'^x^u) du 



We observe that r[s) is an entire function, h{s) is defined and holomorphic for ^s > —b, 
and 

x{s) = Ao{sy^[r{s) + h{s)] (6.5) 

is holomorphic for 9f? s > —b with the exception of finitely many poles. 

Take e > and k > —a. On any strip of the form —b + e < 3? A < k, the 
functions r{s) and h{s) are uniformly bounded. Since k is greater than the real part of 
all singularities of x{s), we can use the inverse formula for the Laplace transform. 



1 



k+ioD 



x{t) = — - / e''x{s)ds. (6.6) 

Choose e > small such that -b + e/2 < k and cr{A) n {s : -b < "Rs < -b + e/2} = 0. 
Note that in the strip —b<^s< —a, the only possible poles of x lie on —b + e/2 < 
3f?s< -a. 

In the strip —6 + e/2 < 3f^s < A;, the functions |r(s)|, |/i(s)| are bounded, and 
||Lo(e*'J)|| < max(l, e^'^"^)||Lo||, hence the operator norm for the inverse of Ao(s) = 
si — Lo(e*'/) satisfies 

l|Ao(s)-i<^ .| , .II, II for |s|>max(l,e-^'OI|i^o||. 

|s|— max(l,e '*^j||Lo|| 

From this estimate and (6.5), we conclude that |e''*a;(s)| — ?■ as 13" s| — )■ oo, uniformly 

in the strip —b + e/2 < 3f^s < A;, and that x{s) is L^-integrable in any straight line 

s = Xo + iy,y eM, for any fixed Xq G [—6 + e/2, k] . We may shift the path of integration 

in (6.6) to the left, and obtain 

x{t) = z(t) + w{t), where 

^(t) = ^i?es(e*-5,A), w{t) = — e'^x{s)ds. (6.7) 

AeA J-b+^~ioo 

From the previous lemma, z{t) is an eigenfunction of (6.2) associated with the set of 
eigenvalues in A. It remains to prove that 

w{t) = 0(e-(^-^)*) at + oo. (6.8) 
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Define u{t) = e^'''''/^^^w{t),v{t) = e(^-')*w(t). We first prove that v E ^^[O, oo). 
Here, Lp[0,oo) = Lp([0, oo); C^),p = 1,2. We have 



M(t) = / e(^+*-^/2)*x(s) ds=— e-'''x{-b + e/2 - is) ds. 

By Plancherel Theorem, ■u(t) G -^^[0, oo), so that v{t) = e~^*/^-u(t) imphes 

||^^IUmo,oo) < ||M||L2[o,oo)||e"'*/lL2[o,oo) < C\\x{-b + e/2 + i-)\\L2(R) < cx), (6.9) 

for some C > 0. Hence, v G L^[0, oo). Define now 

V{t) = Lo{e~^^-'>Vt) = / dr]{e)e~'^^~'^%{t + ^). 
Then, 

/ |\/(t)|rft<max(l,e(^-^)Oll^o||||^||LM-r,oo). (6.10) 

In particular, V G L^[0, oo). We now observe that x{t) = z{t)+w{t), with x{t) a solution 
of (6.1) and z{t) a solution of (6.2). Hence w(t) satisfies (6.1), w'(t) = Lq^Wi) + h(t), 
and we obtain 

v'{t) = {b- e)v{t) + Loie-^'-'^Vt) + e^'-'^'h{t), 

with e(^-^)*/i(t) = 0{e~'^). We conclude that v' G Li[0,oo). Since \v{t)\ < \v{0)\ + 
/o l"^'!-^)! (^^ ^'^^ t > 0, then v is bounded on [0, oo), and (6.8) holds. ■ 

Remark 6.3 For the situation x{t) = 0(e~"*), h{t) = 0(e~^*) (a < 6) as t — )■ oo, denote 
v{t) = {vi{t), . . . ,fjv(t)) as in the above proof. Clearly, one can obtain componentwise 
estimates similar to (6.9) or (6.10). In fact, one concludes that for e > small such 
thata{A)n{s: -b<^s < -b + e/2} = % andt> 0,j = l,...,iV, 



\vj\\l^[o,oo) < (SttVc) \\xj{-b + e/2 + i 



ILV 



and \vj{t)\ < \vj{0)\ + ||-y^-||Li[o,oo) with 

\\Vj\\mo,oo) < C\\xj{-b + e/2 + i-)\\L^R) + ||e^^"'^'/ij||LMo,oo), 

where C = ^ \t^ fi — 'L-^. Similar estimates hold for the case x{t) = 0{e"'*)^h(t) = 
0(e^*) at -oo. 
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